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Summary 

alShin  for  collisiomdominated  plasmas;  Phase  2  sought  to  compare  this  approach  to  published 
Serical  solutions  for  weakly  ionized  plasmas;  and  Phase  3  was  to 

rLent  densitv  Strongly  ionized  plasmas  with  expenments  earned  out  under  other  AfObK 
funding  With  reasonable  success  in  this  final  phase,  the  multifluid  model  was  to  o  provi 
guidance  for  the  design  of  next  generation  arejet  technology. 

The  numerical  multifluid  approach  formulated  in  Phase  1  describ^ 

an  anode  surface  taking  into  account  the  motion  and  interaction  of  electron,  ion,  and  neufral  fields 

boundary.  Ph^  2  a..en.p.s  to  validr^the  dgonto  tev 

not  yet  proven  successful.  Two  different  appro^hes  have  been  “ptor^.  but  ^ve 
demonstrated  numerically  divergent  solutions.  Wrthout  clear  success  m  the  Phase  2  effort, 
comparison  with  experimental  data  in  the  Phase  3  was  postponed. 

Continued  work  with  this  numerical  approach  will  generate  a  physical  correct  and  correlat^ 
aonroach  to  estimating  the  thermal  loading  due  to  the  attachment  of  high  intensity,  high  velocity 
Ts  at tratolphlric  pressures.  The  results  of  this  effort  could  not  only  used  to  design 
arejet  electrodes,  but  also  other  high  current  devices,  such  as  plasma  switches  and  relays. 
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Background 


Low  power  arcjet  thrusters  2  kW)  are  high  specific  impulse,  high  specific  energy  ele^ 
moDubion  devices  that  are  now  used  for  north-south  satellite  stationkeeping.  Arcjet 
Lrformancecanbe  used  to  meet  future  mission  scenarios  if  their  performance 
to  accomplish  other  necessary  functions  such  as  orbit  transfer  or  repositioning.  To  ^ 

adl^ge,  these  additional  on-orbit  r^uirements  -  " 

oneratine  at  significantly  higher  specific  energies  without  a  loss  m  thrust  etticie  y. 
UnfortuLely,  present  data  would  indicate  that  anode  materials  are  now  approachmg  their  design 

limits. 

The  management  of  the  anode  thermal  loading’  is  thus  the  most  critical  area  for  perfonn^ce 

improvement.  An  understanding  of  the  physical  processes  that  govern 

transfer  to  the  anode  is  essential  (particularly  m  the  vicinity  of  the  arcjet  throaO-  Alt  g 

n—  mMeb  are  capable  of  providing  dala  for  critical  design  <^07“. 

depends  on  correctly  modeling  the  underlying  physics  of  arc  aMchment  and  hrat  transfe^ft 

anode-plasma  boundary  -  a  task  most  often  estimated  m  smgle-fluid  approaches  by  electrode- 

wa^models.  The  multifluid  approach  has  potential  to  improve 

providing  separate  descriptions  for  the  electron,  ion,  and  neutral  fields  and  their  exchange  o 
mass,  momentum,  and  energy. 


Objective 


The  objective  of  this  work  was  to  develop  and  verify  analytical  models  wtiich  yield  UB#B  to 
the  key  physical  processes  governing  arc  attachment  and  the  thermal  loadmg  of  anodes  m 

collisional  plasmas. 


Technical  Approach 

The  technical  challenge  of  this  effort  was  to  extend  the  multifluid  analysis  for  weakly  lomzed 
plasmas  of  Meeks  and  Cappelli^  to  significantly  higher  current  densities  and  highly  compressible 

plasmas  typical  of  arcjet  thrusters. 

The  effort  divided  naturally  into  three  phases: 

Phase  1 ;  The  construction  of  a  multifluid  model  of  the  arc  attachment 
process. 
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Phase  2:  Validation  of  the  Phase  1  model  by  direct  comparison  with  the 
recent  work  of  Meeks  and  Cappelh.  Rationalize  differences 
and  refine  the  model  where  necessary. 

Phase  3-  Investigation  of  arc  attachment  phenomena  at  higher  current 

densities  where  significant  coupling  exists  between  the  ^c  and 
bulk  plasma.  This  included  a  companson  of  model  predictions 
with  experimental  data  from  other  AFOSR  research  programs. 

Our  approach  was  to  solve  the  weakly  ionized  stagnation  point 

of  R  eference  2  and  to  directly  compare  the  results  of  the  two  analyses.  The  baseline  test 

rs“"el“siLilartotLofLdinMeek^^ 

toward  the  anode.  The  plasma  constituents  are  Ar  (Class  1),  At  (Class  ),  an  e 

(Class  3).  The  anode  surface  temperature  T^  is  assumed  constant. 


kg  =  Kn^nj  cnfi/  mole-sec 

(surface  recombination  coefficient) 

Figure  1.  Schematic  of  the  Point-Plane  Arc  Attachment  Problem 
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The  sMcification  of  the  incident  current  density  along  this  surface  also  toennines  the  bounty 
^.hton  She  Lgnetic  field  strength  at  the  surface.  In  this  case,  the  neutral  and  ion  cte 
“ImL  were  asstSd  to  obey  the  no-slip  condition.  The  electro^.  howev«,  “"7 
current  and  pass  unhindered  with  a  trajectoiy  normal  to  the  surface.  The  radial  dunension  of  t^ 
attachment  zone,  the  magnitude  of  the  flux  density,  and  a  t^^biMhon^e  mm^^dso^^ 

specified.  The  fiee  stream  conditions  are  taken  to  represent  the  ^  P 
Mvidual  class  temperature,  pressure,  velocity,  and  cheira^  composition 
Unlike  the  approach  in  Reference  2,  the  ion  and  electron  velocities  as  well  as  the  magnehc  field 
are  free  to  change  along  this  surface  as  the  solution  for  current  density  evolves. 

The  center  of  the  arc  column  is  a  line  of  symmetry  for  all  plasma  parameters  with  fte  ^al 

extent  of  the  computational  space  a  distance  R  from  the  centerline.  The  assumption 

Sosen  large  enough,  the  influence  of  the  boundaiy  diminish^  ““  *^2  The  c^oS 

centerline  better  represents  the  one-dimensional  similarity  soluhons  of  ^ 

for  other  parameters  along  the  R-surface,  however,  are  less  obvious.  The  approach  taken  here 

was  to  seUhe  magnetic  field  to  zero  along  this  boundary  to  limit  the  currem  to  the 

space.  The  ion  and  electron  class  boundary  conditions  along  the  surface  were  frea 

continuative  (i.e.,  the  values  of  the  variables  do  not  change  across  the  boun^).  Finally,  the 

neutral  class  pressure  (P)  or  pressure  gradient  parameter  (  defined  as  A  m  Reference  2) 

was  set  to  a  constant  value. 


Results:  Phase  1 


A  multifluid  description  of  the  attachment  process  has  been  developed  to  represent  the  ^ehav 
of  neutrals,  ions,  and  electrons  within  a  plasma.  There  are  no  ^ 

populations  of  each  field.  Each  particle  class  has  its  own  velocity  and  temperature  distribute 
fenerev)  and  is  coupled  to  the  other  classes  by  inter-class  mass,  momentum,  and  ener^  exchange 
processes.  The  derivation  of  these  multifluid  conservation  equations,  auxiliary  state 
relationships,  and  exchange  terms  have  been  discussed  in  previous  work. .  Note  that  m  the 
model  equations  shown  in  Figure  2,  the  magnetic  field  (B)  is  a  primary  variable  and  the  electnc 
field  (E )  and  associated  potential  (<J> )  are  derived  quantities. 


The  Phase  1  work  produced  two  numerical  approaches  to  the  solution  of  these  equations  and  the 
construction  of  the  associated  algorithms  (i.e..  Algorithm  1  and  Algorithm  2).  The  two 
approaches  differed  fundamentally  in  the  construction  of  the  solver  for  O. 


Algorithm  1  4>- Solver 

Construction  of  the  Algorithm  1  d>-solver  was  completed  early  in  the  program  and  initially 
validated  by  comparison  of  its  solution  with  a  number  of  readily  available  static  test  cases. 
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Incorporation  of  the  ^.-solver  with  the  parent  multifluid  was  straightforward  since  it  was 
constructed  to  operate  on  the  same  mesh  as  that  used  for  the  plasma. 

andVihon.  followed  by  ~  “r^p^tlfl^  JctZ  "  and  S,. 

fields.  Solving  Equation  6  for  B  wift  the  electron  velocity.u.,  ^sumed  ^ 

solution  for  the  total  current  density  J  (Equation  5),  Subtrac^  the  ton  cuirent  J  le^ 
directly  to  the  electron  current  and  an  updated  value  for  u.  With  the  updated  eshmates  fo 

tlociL,  and  additional  sources  of  energy  transport  calculat^  ‘^TerenS'^iSS’ 
Ohmic  hatting,  etc.),  the  energy  transport  (Equation  3)  is  solved  to  the  ” 

At  this  Doint  the  exchange  of  mass,  momentum,  and  energy  between  Ae  difterent  cias 

accomplished  using  appropriate  collision  cross  sections  and  ^  ^ 

Ugra.^  state  variables  are  mapped  firom  the  Lagran^  space  back  onto  the  Eulenan  gnd 

This  cycle  is  repeated  until  a  steady  state  solution  is  obtained. 

Static  plasma  calculations  (i.e.,  all  initial  field  velocities  set  to  zero)  demonstrated  the  anticipated 
null-solutions  without  spurious  numerical  errors.  Differences  from  eye  e-to-cyc  e  were  o 
order  of  the  anticipated  numerical  round  off.  However,  for  plasmas  with  non-zero 
numerical  solution  slowly  diverged,  demonstrating  a  continually  incr^mg 
or  non-neutral  plasma  condition.  This  scenano,  of  course,  violates  the  charp  neutrality 
assumption  upon  which  the  algorithm  is  based  and  led  to  re-evaluation  and  eventual  red^ign  of 
the  ni^erical^proach.  The  refined  methodology.  Algorithm  2,  is  discussed  m  the  following 

section. 

An ev^1^n"^fthr Algorithm  1  solution  suggested  the  need  f  f 

solution  method.  The  principle  problem  exposed  by  the  approach  taken  m  Algonthm  1 

the  prediction  for  the  electron  velocity  depended  on  the  ion  velocity  estimate  which,  m  turn,  wa 

derived  from  the  total  current  density.  Note  that  the  current  density  depends  on  B  *®  ^ 
estimate  for  ue .  An  error  in  any  of  one  of  these  steps  leads  to  further  propagation  of  errors  and 

decreasing  solution  fidelity. 

An  alternate  solution  (Algorithm  2)  method  suggested  by  Brackbill^  makes  use  of  the  equivalent 
quasi-neutrality  condition  VJ=  0  where  J  is  the  total  current,  to  provide  an  indirect 
calculation  of  the  electron  current.  Followmg  Brackbill,  the  mertial  terms  m  the  electro 
momentum  equation  (Equation  2)  are  neglected  ( P.due/dt «  0 )  and  an  equation  for  the  electron 
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current  density  is  then  derived  by  multiplying  Equation  2  by  q/m,  where  q,  is  the  electron 
charge  and  rck  is  the  electron  mass; 


T  _  =.^E  +  5. 

‘  Ve  n. 

«  is  the  electron  plasma  frequency  squared,  q,  =  is  the  electron 

m. 

qK  ...  . 

cyclotron  frequency,  i/*  =  2  resistivity,  and 


.1^ 


where  0)%- 


cnte 


_  1 
^‘-7, 


— (-V  +  V  •  re)  +  1-2  KteJt 

me 


Following  Brackbill’s  derivation  and  writing  expressions  for 

X  Q,  and  J,  •  a  to  solve  for  j, ,  leads  to  an  expression  for  the  electron  current  density. 


J  =  - L_—  .  \(Ope\  E  +  E  X  +  (E  •  Q  ^Qe]  +  [Se  +  5e^  Qe  +  (^e  * 

(l  +  de)  ^ 

where  and  Q,  =  ^  •  Using  V  •  J  =  0  =  V  •  (j,+  J  J  and  assuming  that  j,  is 

Ve  Ve  .  -  U  A  ■  A. 

known  from  the  previous  ion  momentum  equation,  an  equation  for  E  can  be  denved. 


V(e-E)  =  -V*l-Ij,. 

'  '  t^e 

This  equation  can  be  solved  for4.assummgthat-i-f.VxE-  0  which  in  turn  leads 

to  the  approximation  E  =  -  Vd> . 

The  logical  flow  of  Algorithm  2  then  follows: 

1 .  Solve  for  the  rate  of  change  of  plasma  species,  P, 

2.  Solve  the  momentum  equation  for  both  ions  and  neutrals 

3.  Evaluate  the  ion  current,  j, 

4.  Solve  for  the  potential  field  and  the  electric  field 

5 .  Solve  for  the  electron  velocity 

6.  Solve  for  the  magnetic  field  using  Ampere’s  Law 

7.  Advance  the  neutral,  ion,  and  electron  energy  equations 
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8.  Return  to  the  top  of  the  cycle  and  begin  again. 

Although  this  scheme  proved  more  difficult  to  implement  than  ^  “f  ‘j 

solutions  obtained  for  flowing  plasmas  with  simple  boundary  conditions  were  stable 
remained  quasi-neutral  regardless  of  calculation  time. 

Comparison  of  the  Algorithm  2  results  with  the  weakly  ionized  solutions  of  R®[erenc®  2, 
however,  again  proved  unsuccessful,  largely  due  to  difficulties  in  implenientmg  the  more 
complex  boundary  conditions.  Small  errors  introduced  by  the  numencal  approach  at  the 
boundaries  propagated  into  the  computational  space,  eventually  corrupting  the  solut  . 
A  complete  investigation  of  this  problem  and  possible  solutions  were  not  completed 
before  the  end  of  the  present  contract  period. 


Conclnsions 

This  work  has  exposed  shortcomings  with  standard  multifluid  approaches.  It  has  also 
investigated  an  promising  approach  using  a  method  outlined  by  Brackbill  to  solve  t  e 
multifluid  equations.  Further  work  is  required  to  understand  and  implement  boundary 
conditions  in  this  new  scheme.  The  completion  of  this  task,  however  is  ^ential  if  ^ 
are  to  compare  the  more  general  multifluid  approach  with  the  work  of  Reference  2  or  to 
gain  further  insight  into  the  physics  of  arc  attachment. 
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Directions  for  Arcjet  Technology  Development 


G.  W.  Butler*  and  R.  J.  Cassadyt 
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Guidelines  for  improvements  in  arcjet  performance  based  on  a  simplified  analysis  of  advanced  mission 
requirements  are  presented.  The  results  of  this  analysis  clearly  show  that  to  remain  competitive  with 
next-generation  electric  propukion  devices,  significant  increases  in  the  specific  impulse  of  present  arcjet 
technology  must  be  achieved.  Of  the  various  design  options  available,  concepts  that  improve  the  thrust 
efficiency  offer  the  greatest  performance  leverage,  and  may,  in  fact,  represent  the  only  approaches  capable 
of  meeting  the  advanced  mission  goals.  If  they  are  to  be  successful,  such  concepts  will  improve  the 
management  of  energy  and  mass  transport  within  the  arcjet.  One  method  to  understand  and  control  the 
internal  plasmadynamics  of  the  arcjet  to  improve  its  performance  is  outlined. 


Nomenclature 


H* 

/sp 

m 

NjH4 

P 


Pc 

Pn 

Palm 

TIP 

W-100 

W-2%Th02 

W-4%Re-0.3%HfC 

AV 

V 

Vpcu 

cr. 


=  excited  state  hydrogen  atom 
=  specific  impulse,  s 
=  propellant  mass  flow  rate,  mg/s 
=  monopropellant  hydrazine 
=  power  delivered  to  the  propulsion 
system,  W 

=  chamber  pressure,  MPa 
=  power  delivered  to  the  propellant 
(plasma),  W 
=  specific  energy,  J/kg 
=  thrust-to-power  ratio,  N/W 
=  pure  tungsten  metal 
=  thoriated  tungsten 
=  tungsten -rhenium- hafnium  carbide 
alloy 

=  velocity  increment,  m/s 
=  thrust  efficiency,  % 

=  power  conditioning  unit  conversion 
efficiency,  % 

=  electrical  conductivity 


Introduction 

Arcjet  thrusters  are  high  specific  impulse,  high  specific 
energy,  electric  propulsion  (EP)  devices  now  in  use  for 
North -South  (N-S)  satellite  stationkeeping  of  geosynchro¬ 
nous  satellites.  At  present,  low-power  arcjets  operate  on  hy¬ 
drazine  decomposition  products  and  afford  an  /,p  in  the  500— 
600-s  range  with  tj  of  approximately  30-35%.  Figure  1 
presents  a  schematic  of  the  physics  at  work  within  a  conven¬ 
tional  arcjet  thruster.  The  discharge  is  a  wall-  and  vortex-sta¬ 
bilized  arc.  The  arc  emanates  from  the  cathode,  passes  through 
the  constrictor,  and  attaches  in  a  diffuse  (nondestructive)  man¬ 
ner  along  the  expanding  portion  of  the  nozzle  (anode).  The 
propellant  is  injected  upstream  of  the  cathode  tip  with  an  azi¬ 
muthal  component  of  velocity.  The  resultant  swirl  promotes 
arc  stability  and  acts  to  reduce  heat  transfer  to  the  upstream 
anode  surfaces.  The  propellant  merges  with  the  arc  as  it  passes 
from  the  plenum  through  the  constrictor,  into  the  expanding 
portion  of  the  nozzle,  and  finally  to  ambient  space.  The  merg¬ 
ing  of  the  arc  with  the  propellant  gas  is  a  complex  and  con¬ 
tinuous  process.  No  distinct  boundary  between  arc  and  nonarc 
fluid  can  be  drawn. 
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The  shaded  zones  of  Fig.  1  identify  regions  of  the  discharge 
in  which  different  physical  processes  may  be  considered  dom¬ 
inant.  The  principal  activity  in  zone  1  is  the  dissipation  of 
electrical  energy  via  collisions  of  electrons  with  the  neutral 
and  ionized  propellant.  Some  of  the  ions  produced  in  this  re¬ 
gion  migrate  to  the  cathode  with  sufficient  energy  to  maintain 
the  discharge.  The  bulk  of  the  current  is  made  up  of  electrons 
that  travel  through  the  propellant  under  the  influence  of  the 
electric  field  to  the  anode.  They  are  continually  decelerated  by 
collision  with  their  heavier  neighbors  and  continually  accel¬ 
erated  by  the  imposed  electric  field.  In  zone  2,  this  process 
results  in  a  substantial  increase  in  the  bulk  temperature  of  the 
plasma.  This  situation  changes  dramatically  as  the  plasma 
leaves  the  constrictor  and  enters  zone  3,  which  is  characterized 
by  decreasing  field  strengths  and  rapid  gasdynamic  expansion 
of  the  plasma.  Thermal  energy  exchange  between  the  electrons 
and  other  species  thus  falls  dramatically  so  that  this  region  is 
dominated  by  the  evolution  of  energy  from  recombination  re¬ 
actions.  In  zone  4,  the  conversion  of  this  internal  energy  to 
kinetic  energy  occurs  through  gasdynamic  expansion  pro¬ 
cesses.  In  this  zone,  viscous  dissipation  and  heat  transfer  at 
the  anode  surface  reduce  the  benefits  of  the  expansion  pr^ess. 

Typical  low-power  arcjets  have  conical  converging -diverg¬ 
ing  nozzles  with  constrictor  diameters  on  the  order  of  0.05  cm, 
nominal  expansion  angles  of  20  deg,  and  exit  diameters  of  0.35 
cm.  Moving  up  in  power  to  30  kW,  representative  constrictor 
and  exit  dimensions  increase  to  0.5  and  2.5  cm,  respectively. 
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CONVERSION  OF  ELECTRICAL  TO  CHEMICAL  POTENTIAL 
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CONVERSION  OF  CHEMICAL  POTENTIAL  TO  INTERNAL  ENERGY 
CONVERSION  OF  INTERNAL  ENERGY  TO  KINETIC  ENERGY 

Fig.  1  Key  features  of  the  arcjet  physics. 
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The  physical  characteristics  of  the  arcjet  fiowfield  vary  from 
a  nearly  fully  ionized  plasma  at  temperatures  in  excess  of 
20,000  K  near  the  cathode  tip  to  a  relatively  cold  plasma 
(1000-2000  K)  at  the  anode  surface.  Moreover,  velocities  vary 
from  approximately  10  km/s  on  centerline  to  zero  at  the  wall. 
Considering  the  magnitude  of  these  parameters  and  small  ge¬ 
ometric  length  scales,  it  is  clear  that  the  arcjet’s  performance 
depends  largely  on  nonequilibrium  transport  processes. 

Future  mission  scenarios  discussed  in  the  following  section 
reflect  the  need  for  greater  control  of  on-orbit  assets.  The  sta¬ 
tionkeeping  roles  of  present  low  power  arcjets  must  expand  to 
include  orbit  maneuvering  and  some  limited  number  of  repo¬ 
sitioning  events.  For  the  arcjet,  this  will  require  even  higher 
/,p,  which,  in  turn,  means  driving  their  operational  envelope  to 
higher  specific  energies  (Pn/'fi)  while  maintaining  or  improv¬ 
ing  thrust  efficiency. 

The  present-day  technology,  however,  is  limited  by  the  an¬ 
ode’s  ability  to  withstand  the  high  thermal  fluxes  that  are  char¬ 
acteristic  of  operation  at  high  Palni  (e.g.,  higher  power  dissi¬ 
pation  and/or  lower  mass  flow  rates).  In  fact,  many  prototype 
thrusters  fail  in  extended  life  tests  because  of  the  intolerable 
heat  loads  on  the  anode.  A  significant  performance  improve¬ 
ment  without  loss  of  useful  life  must  therefore  be  achieved 
through  some  combination  of  higher  higher  /,p,  or  higher 
Tj.  The  motion  of  the  plasma,  the  large  gradients  in  chemical 
composition  and  fluid  properties,  and  the  relatively  short  res¬ 
idence  times  make  this  a  formidable  design  problem. 

Advanced  Mission  Requirements 

The  mission  applications  for  on-orbit  propulsion  can  be  gen¬ 
erally  collected  into  three  categories:  1)  stationkeeping,  2)  ma¬ 
neuvering,  and  3)  transfer.  These  functions  were  referred  to  in 
a  recent  study  even  more  simply  as  1)  hold,  2)  move,  and  3) 
lift.'  We  will  describe  each  of  these  functions  and  the  impli¬ 
cations  of  improvements  in  arcjet  performance.  Arcjet  perfor¬ 
mance  for  both  a  baseline  and  an  advanced  technology  level 
will  be  compared  with  other  EP  thruster  systems.  Table  1  lists 
the  parameters  assumed  for  each  propulsion  system  separated 
into  categories  of  hold  and  move  and  lift.'^  ^  The  system  power 
is  defined  as  the  power  to  the  power  conditioning  unit  (PCU) 
and  the  life  and  total  impulse  numbers  are  the  assumed  capa¬ 
bilities  of  a  single  thruster. 

Hold 

Hold  maneuvers  maintain  a  spacecraft  in  a  given  orbit.  Ex¬ 
amples  include  (N-S)  and  East-West  (E-W)  stationkeeping 
of  geosynchronous  (GEO)  communications  satellites  and  drag 
makeup  for  low  earth  orbit  (LEO)  satellites.  Thrust  levels  are 
not  a  primary  consideration  because  the  perturbing  force  is  of 
low  magnitude.  Higher  Ap  is  the  more  important  performance 
parameter  since  the  fuel  mass  required  for  the  mission  must 
be  minimized.  Electric  propulsion  is  well  suited  to  this  require¬ 


ment.  In  fact,  electric  propulsion  devices  such  as  pulsed 
plasma  thrusters  (PPTs)  and  electrothermal  hydrazine  thrusters 
(EHTs)  have  been  routinely  used  since  the  early  1980s  for  drag 
makeup  and  N-S  stationkeeping,^’  respectively.  In  early  1994, 
arcjets  were  used  for  the  first  time  aboard  the  Telstar  401  sat¬ 
ellite  for  N-S  stationkeeping. 

The  EP  systems  described  in  Table  I  were  evaluated  for  the 
N-S  stationkeeping  mission  of  a  GEO  satellite  with  an  initial 
mass  of  1500  kg.  Satellite  lifetime  was  varied  from  10  to  15 
years.  Power  delivered  to  the  propulsion  systems  was  assumed 
to  be  limited  to  4  kW.  Thruster  concepts  were  assumed  to  use 
the  demonstrated  power  level  for  the  particular  device  listed 
in  Table  1  and  the  number  of  thrusters  was  scaled  to  come  as 
close  as  possible  to  4  kW.  Using  parametric  analysis,  the  total 
propulsion  system  mass  was  evaluated  as  a  function  of  pay- 
load,  The  EHT  (300-s  MR-502)  system  mass  is  highest  across 
the  range  of  AV  required  for  a  10-  15-year  mission.  The  base¬ 
line  arcjet  (502-s  MR-508)  falls  in  an  intermediate  range.  An 
improved  arcjet  (650  s)  is  much  more  competitive  on  a  system 
wet  mass  basis  with  the  30-cm  ion  engine  and  the  SPT-100 
Hall  thruster. 

The  second  consideration  in  the  application  trades  is  the 
time  required  to  perform  the  stationkeeping  maneuver.  This  is 
inversely  proportional  to  the  thrust  level  of  the  EP  device.  For 
a  fixed  power  level  delivered  to  the  propulsion  system,  this 
trade  favors  a  device  with  a  high  TIP,  as  shown  in  Fig.  2.  It 
is  also  important  to  maintain  or  improve  thrust  efficiency  as 
shown  by 


(1) 

2Pft  2niPfi  2  \P n/  2i7pcu  \P) 

Reduced  thrust  efficiency  requires  longer  or  more  frequent 
thruster  bums  to  perform  the  hold  mission,  complicating  bat¬ 
tery  and  power  distribution  system  designs.  In  general,  longer 
bums  are  less  desirable  to  users  because  mission  operations 
(ground  control)  become  more  complex.  Another  consideration 
is  the  life  testing  that  drives  the  cost  to  qualify  the  system. 
Typical  users  require  the  demonstration  of  lifetime  1.5  times 
the  mission  life.  Obviously,  systems  that  require  longer  thruster 
life  will  cost  more  to  qualify. 

Figure  3  shows  approximate  operational  envelopes  for  cur¬ 
rent  and  advanced  arcjet  technology.  Increasing  /,p  in  the 
l.5-2.0-kW  power  range  is  critical  to  hold  mission  applica¬ 
tions.  This  will  make  arcjet  systems  mass  competitive  with  the 
advanced  Hall  and  ion  thmster  systems  for  payloads  up  to 
1500  kg  and  mission  lifetimes  of  up  to  15  years.  Because  of 
the  user  concern  over  operational  issues,  it  is  also  important 
to  maintain  rj  as  Ap  is  increased  to  minimize  the  frequency  and 
duration  of  thmster  firings. 


Table  1  EP  systems  summary 


Name 

System 

power, 

W 

Thrust, 

N 

s 

Life, 

h 

Inert 

mass, 

kg 

Tank 

fraction 

Total  impulse 
capability, 
N-s 

Hold  and  move  propulsion  systems 

MR-502  EHT 

750 

3.90E-01 

302 

370 

1.85 

0.04 

5.19E+05 

MR-508  arcjet 

1800 

2.25E-01 

502 

870 

7.50 

0,04 

7.05E+05 

650-s  arcjet 

2000 

2.01  E-01 

650 

900 

7.50 

0.04 

6.51E+05 

XIPS-13 

500 

1.78E-02 

2585 

6602 

16,98 

0.20 

4.23E+05 

NASA  30-cm  ion 

1822 

7.30E-02 

2814 

8000 

23.50 

0.20 

2.10E+06 

SPT-100 

1500 

7.01E-02 

1650 

4000 

16.30 

0.20 

1 ,01E+06 

Lift 

propulsion 

systems 

H2  arcjet 

10,000 

6.0E-01 

1200 

1500 

20 

0.10 

3.24E+06 

Ion 

2500 

1.3E-01 

3200 

8000 

30 

0.15 

3.74E+06 

SPT/Hall 

10,000 

8.2E-01 

1500 

4000 

25 

0.15 

1.18E+07 
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Fig.  2  Comparison  of  thruster  bum  time  for  hold  maneuvers. 


Fig.  3 


800j 

750  4 


0.20  0,30  0.40  0.50 

Thrust  Efficiency 


1  Target  Arcjet  Performance 
Present  Arcjet  Performance 
Advanced  arcjet  performance  for  hold  and  move  missions. 


Move 

On-orbit  maneuvering  of  satellites  is  of  increasing  interest 
to  users.  This  category  encompasses  orbit  adjustment  (e.g.,  cir¬ 
cularization),  plane  changes,  and  E-W  reposiUoning,  among 
others.  In  the  past,  satellite  assets  have  been  launched  and 
inserted  into  orbits  that  were  more  or  less  fixed  for  the  life  of 
the  satellite.  New  demands  for  lower  cost  systems  are  driving 
users  toward  a  different  strategy,  employing  fewer  satellites 
that  are  capable  of  moving  as  needs  change.  For  rapid  repo¬ 
sitioning  (i.e.,  30  days  or  less),  the  AV  required  for  one  180- 
deg  repositioning  maneuver  is  approximately  the  same  as  1 
year  of  N-S  stationkeeping.  The  use  of  EP  systems  thus  en- 
ables  such  missions. 

A  typical  mission  in  this  category  is  E-W  repositioning  of 
a  GEO  satellite.  Figure  4  shows  the  results  of  a  comparison 
of  total  propulsion  system  mass  between  candidate  EP  systems 
for  an  E-W  maneuver  halfway  around  the  equator.  To  enable 
a  direct  comparison,  a  maneuver  time  of  21  days  was  selected. 
Power  available  for  the  EP  systems  was  assumed  to  be  limited 
to  2.0  kW  because  a  typical  military  satellite  has  less  onboard 
power  than  similar  communications  satellites.  Actual  power 
for  each  system  was  determined  by  their  demonstrated  capa¬ 
bilities  as  shown  in  Table  1 . 

The  results  show  the  effect  of  the  lower  inert  mass  of  the 
arcjet  system  trading  against  the  higher  /^p  of  the  ion  and  Hall 
thrusters.  The  EHT  and  XIPS-13  thrusters  are  not  competitive 
on  a  system  mass  basis  for  this  mission  and  are  not  discussed. 
The  SPT-100  Hall  thruster  system  weighs  less  than  the  MR- 
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Fig.  4  Comparison  of  propulsion  system  mass  for  move  missions. 


508  arcjet  over  the  entire  range.  The  30-cm  ion  engine  trades 
poorly  with  the  MR-508  arcjet  below  five  repositioning  ma¬ 
neuvers,  but  has  a  mass  advantage  if  more  than  six  maneuvers 
are  desired.  The  advanced  arcjet  design  (650  s)  surpasses  the 
ion  engine  and  is  competitive  with  the  Hall  thruster  over  the 
entire  range.  While  not  shown  here,  these  trade  studies  have 
also  revealed  the  importance  of  high  arcjet  thrust  efficiency  for  ^ 
time-critical  maneuvers.  Although  this  is  of  less  consequence 
for  the  assumed  mission  duration,  it  is  still  of  concern,  espe¬ 
cially  if  the  desired  maneuver  time  is  much  less  than  the  as¬ 
sumed  21  days. 

The  approximate  move  mission  performance  ranges  for  an 
N2H4  arcjet  are  the  same  as  the  hold  maneuver  (Fig.  3).  The 
main  reason  to  select  an  arcjet  for  on-orbit  maneuvering  is 
once  again  increased  Ap.  The  thrust  efficiency  must  be  main¬ 
tained  for  the  arcjet  to  retain  its  bum-time  advantage,  partic¬ 
ularly  in  situations  where  total  maneuver  durations  of  10-14 
days  are  desired. 

Lift 

The  final  mission  category  examined  is  orbit  transfer  or  lift. 
Mission  planners  have  long  evaluated  the  concept  of  an  elec¬ 
tric  orbit  transfer  vehicle  (EOTV).*"'®  EOTVs  have  the  major 
advantage  of  large  propellant  mass  savings  relative  to  chemical 
upper  stages.  The  major  drawback  to  EOTV  systems  is  the 
long  trip  times  required  to  complete  the  transfers.  Preventing 
the  trip  times  from  becoming  excessively  long  requires  that 
power  levels  for  the  EOTV  increase  to  the  30-50-kW  range. 
High-power  levels  such  as  these  require  large  (200-300  m  ) 
thin  blanket  deployable  solar  arrays,  which  currently  can  cost 
up  to  two  million  dollars  per  kilowatt.  Minimizing  the  power 
requirement  is  therefore  desirable  to  maintain  a  reasonable  cost 
and  size  for  the  solar  arrays. 

A  typical  lift  mission  example  is  the  deployment  of  a  de¬ 
fense  support  program  (DSP)  satellite  (2500  kg)  to  GEO  from 
a  LEO  parking  orbit.  Figure  5  shows  the  results  of  a  trade 
between  EP  thruster  systems  for  this  mission.  The  thruster  per¬ 
formance  parameters  are  listed  under  the  Lift  Propulsion  Sys¬ 
tems  heading  in  Table  1.  EOTV  power  level  was  assumed  to 
be  30  kW  for  all  concepts.  The  EP  systems  all  save  consid¬ 
erable  propellant  mass  over  the  chemical  system.  Trip  times 
shown  in  parentheses  indicate  an  advantage  for  the  arcjet  be¬ 
cause  of  its  higher  TtP.  The  hydrogen  (H2)  arcjet  must  achieve 
an  Ap  of  1300  s  and  a  thrust  efficiency  of  40%  to  be  mass 
competitive  and  maintain  its  trip  time  advantage.  This  repre¬ 
sents  the  minimum  acceptable  performance  and,  as  shown  in 
Fig,  6,  is  not  within  the  range  of  current  performance  capa¬ 
bility.  H2  arcjet  technology  for  EOTV  missions  is  the  area  with 
the  largest  gap  between  current  capability  and  required  per¬ 
formance. 
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Fig.  5  Payload  mass  comparison  for  DSP-class  satellite  lift  maneuver. 
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Fig.  6  Hydrogen  arcjet  performance  target  for  lift  missions. 


Another  example  of  a  lift  mission  is  orbit  acquisition  for  a 
LEO  constellation.  Examples  of  such  constellations  are  the 
Iridium  network  and  the  Teledesic  system  now  in  concept  def¬ 
inition."  These  constellations  involve  large  numbers  of  satel¬ 
lites  in  many  orbit  planes.  The  initial  insertion  into  a  parking 
orbit  is  accomplished  by  putting  4-10  satellites  on  a  single 
launch  vehicle.  From  the  parking  orbit,  each  satellite  must  then 
be  raised  to  its  final  altitude  and  moved  to  its  operational  po¬ 
sition.  Assumed  performance  parameters  for  the  thruster  sys¬ 
tems  are  the  same  as  for  hold  and  move  missions.  System 
power  was  assumed  to  be  5.0  kW. 

Figure  7a  and  7b  show  the  results  of  a  comparison  of  EP 
devices  for  a  representative  LEO  acquisition  mission,  such  as 
the  deployment  of  a  700-kg  satellite  to  an  orbital  altitude  be¬ 
tween  500-1000  km.  The  results  show  that  the  arcjet  Ap  must 
be  improved  to  be  mass  competitive  with  the  Hall  thruster  and 
ion  engine.  The  arcjet  TIP  again  provides  a  bum  time  advan¬ 
tage  over  all  the  other  EP  concepts.  Figure  8  shows  the  de¬ 
velopment  emphasis  for  this  class  of  arcjet  thruster  is  toward 
improved  while  maintaining  or  slightly  increasing  t]. 


Arcjet  Development  Strategies 
The  improvements  required  to  meet  the  next-generation  per¬ 
formance  goals  can  be  realized  only  if  new  designs  are  able 
to  achieve  some  measure  of  control  over  the  key  physical  pro¬ 
cesses  described  earlier  and  shown  in  Fig.  1.  Figure  9  illus¬ 
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Fig.  7  Summary  of  electric  thruster  system  trades  for  LEO  orbi 
acquisition:  a)  dependence  of  system  mass  and  b)  burn  time  oi 
AV  requirements. 


trates  a  generic  arcjet  performance  map  that  will  serve  as  back 
ground  for  this  discussion.  Each  curve  represents  the  variatioi 
of  1)  and  Ap  with  m  for  one  value  of  Pn-  Pn/w  increases  anc 
m  decreases  along  each  curve  in  the  direction  of  increasing  Aj 
Note  that  a  maximum  value  for  h  is  shown  for  each  curve 
Recent  Hz  arcjet  experiments  have  shown  that  Ap  is  strong  1; 
influenced  by  frozen-flow  losses  below  this  maximum,  and  b; 
thermal  losses  above  it.*^  As  Pn/m  increases,  the  time  scale 
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Fig.  8  Arcjet  performance  target  for  LEO  orbit  acquisition. 


Fig.  9  Generic  arcjet  performance  map  showing  three  paths  to 
achieve  higher 


associated  with  diffusion  and  thermal  conduction  become  sig¬ 
nificant  compared  to  the  plasma  transit  time.  Operation  at  high 
Pnfrh-high  is  thus  limited  largely  by  thermal  rather  than 
frozen-flow  losses.  Above  some  critical  value  of  Pn/'w,  the 
energy  flux  to  the  anode  severely  degrades  thruster  efficiency 
and  reduces  thruster  life.  It  follows  that  a  critical  value  for 
PJrh  can  be  defined,  at  least  qualitatively,  in  terms  of  the 
maximum  steady-state  energy  flux  that  the  anode  can  tolerate. 

Consider  the  Ap  value  indicated  by  the  solid  line  in  Fig.  9 
as  an  advanced  hold  or  move  mission  goal  and  curve  Pi  to  be 
representative  of  present  thruster  performance  at  constant 
power.  Three  design  options  exist  to  provide  higher  /,p  and 
greater  mission  flexibility: 

1)  Maintain  Pn  and  decrease  m  following  curve  P,  to  reach 
the  Ap  goal.  This  is  the  short-term  arcjet  development  strategy. 
Higher  performance  is  being  pursued  by  simply  reducing  m  at 
constant  Pn-  The  success  of  this  approach  depends  upon  the 
development  of  anode  materials  that  are  capable  of  withstand¬ 
ing  significantly  higher  energy  fluxes. 

2)  Increase  P^  to  the  Pj  value  along  a  line  of  constant  Pn/ 
m  and  then  decrease  m  to  reach  the  /^p  goal.  This  option  re¬ 
quires  operating  the  arcjet  at  higher  Pn  and  higher  m.  Note 
that  increases  in  thruster  power  would  still  be  accompanied  by 
a  reduction  in  m  to  attain  the  design  goal. 

3)  Maintain  Pn  (P3  =  Pi),  increase  17  via  hardware  improve¬ 
ments,  and  decrease  rh  to  reach  the  /,p  goal.  This  is  the  most 
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Fig.  10  State-of-the-art  and  advanced  arcjet  performance  re¬ 
quirements  for  hold  and  move  missions. 

difficult  of  the  three  strategies  since  certain  key  physical  pro¬ 
cesses  must  be  controlled  to  produce  higher  17s. 

Figure  10  illustrates  the  state-of-the-art  performance  levels 
of  low-power  N2H4  thrusters  for  hold  and  move  missions.  The 
graph  shows  the  variation  of  /sp  with  77  for  the  NASA  Ad¬ 
vanced  Technology  Development  thruster  operating  at  1 .8  kW 
and  the  Olin  Aerospace  Model  MR-509  thruster  operating  at 
1.6  kW.  The  shaded  box  represents  the  advanced  mission  re¬ 
quirements  outlined  in  the  previous  section.  Extrapolation  of 
these  data  to  higher  /,p  makes  it  clear  that  simply  increasing 
Palm  (option  1)  will  not  suffice.  It  is  also  unlikely  that  onboard 
power  available  for  advanced  low-power  thrusters  will  signif¬ 
icantly  increase.  This  means  that  simply  increasing  the  bus 
power  (option  2)  will  not  provide  any  appreciable  benefit.  Op¬ 
tion  2  is  still  viable  for  EOTV  missions  since  power  level, 
lifetime,  and  configuration  (e.g.,  the  use  of  two  10-kW  or  four 
5-kW  thrusters)  are  somewhat  flexible.  Trip  time  constraints 
and  the  cost  of  solar  arrays,  however,  will  drive  the  EOTV 
configuration  toward  the  highest  achievable  TIP.  From  these 
considerations,  it  is  evident  that  improving  thruster  designs  to 
increase  operating  efficiency  (option  3)  is  the  best  strategy  to 
achieve  the  next-generation  performance  and  lifetime  goals. 
This  approach  also  affords  the  greatest  technology  leverage 
since  high- 17  design  concepts  can  be  applied  to  a  broad  range 
of  arcjet  power  classes  without  affecting  their  TIP  advantage. 

Management  of  Energy  Transport 
Within  the  Anode 

The  choice  of  arcjet  anode  materials  represents  the  true  start¬ 
ing  point  for  discussions  of  performance  improvement.  Any 
improvement  in  the  efficiency  of  the  basic  plasma  processes 
will  directly  add  to  benefits  derived  from  new  electrode  ma¬ 
terials  or  from  novel  schemes  to  process  the  electrode  energy 
flux.  Performance  requirements  force  the  anode  material  to  op¬ 
erate  at  the  highest  Pn/^  can  tolerate  over  the  thruster’s 
useful  life.  The  data  of  Fig.  1 1  show  the  variation  of  {PJtn)^^^ 
(a  measure  of  the  inverse  of  constrictor  diameter)  with  time 
for  life  tests  of  a  conventional  low-power  N2H4  arcjet  operating 
at  1.8  kW  and  mission-averaged  /^p  of  600  s.  The  symbols 
represent  the  behavior  of  the  same  anode  geometry  for  three 
different  materials:  W-100,  W-“2%Th02,  and  W-4%Re- 
0.3%HfC,  The  solid  line  indicates  an  acceptable  rate  of  in¬ 
crease  for  this  parameter  derived  from  the  flight  qualification 
data  of  MR-508,  operating  at  1.6  kW  and  approximately  500- 
s  /jp. 

The  underlying  problem  identified  by  these  data  is  not  that 
the  anode  fails,  but  rather  that  TIP  falls  off  as  the  constrictor 
diameter  decreases,  and  eventually,  the  thruster  no  longer  sat¬ 
isfies  the  mission  requirements.  For  example,  increasing  the 
power  from  1.6  to  1.8  kW  produces  large  increases  in  (P^/m)*'^ 
for  both  pure  W-100  and  W-2%Th02,  and  greatly  reduces 
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Fig.  11  Arcjet  chamber  pressure  variation  caused  by  anode  ther¬ 
mal  effects. 

the  probability  of  flight  qualification.  The  constrictor  closure 
rate  of  the  W-4%Re-0.3%HfC  anode,  however,  falls  below 
and  even  improves  upon  that  of  present  flight  qualified  de¬ 
signs.  This  demonstrates  that  improvements  in  electrode  ther¬ 
mal  properties  can  permit  operation  at  higher  /^p  or  Pnim  with 
considerable  gains  in  performance  and  life. 

In  addition  to  the  development  of  new  anode  materials,  other 
ways  to  efficiently  process  higher  thermal  fluxes  must  also  be 
developed.  An  obvious  choice  is  to  simply  cool  critical  regions 
of  the  anode  with  the  incoming  propellant.  For  example,  con¬ 
sider  the  benefit  of  providing  regenerative  passages  in  a  10- 
kW  Hz  arcjet  thruster.  A  gain  of  nearly  five  percentage  points 
in  7)  at  10  kW  with  ranging  from  100  to  175  MJ/kg  has 
been  reported  in  recent  work.‘^  The  maximum  measured  effi¬ 
ciency  was  approximately  40%  at  Fn/'h  =150  MJ/kg.  A  ther¬ 
mal  analysis  using  measured  temperatures  indicated  that  the 
regenerative  efficiency  achieved  in  this  test  series  was  approx¬ 
imately  94%. 

It  is  our  opinion  that  improvements  in  anode  materials  or 
energy  processing  capability  will  not  by  themselves  permit  op¬ 
eration  at  the  performance  levels  desired  for  advanced  appli¬ 
cations.  Also,  a  question  remains  whether  regenerative  con¬ 
cepts  can  be  successfully  applied  to  low-power  geometries, 
primarily  because  of  their  small  size  and  high  inlet  propellant 
temperatures. 

Management  of  Plasma  Energy  and  Mass  Transport 

Estimates  for  the  distribution  of  electrical  energy  among  the 
various  plasma  processes  described  in  Fig.  1  are  shown  in  Fig. 
12.  The  charts  are  drawn  for  Hz  arcjets  operating  at  a  nominal 
Pn  of  1-5  and  10  kW.  The  values  for  each  category  are  taken 
from  numerical  performance  estimates. Note  that  the  thrust 
energy  does  not  change  significantly  even  when  Pn  increases 
by  a  order  of  magnitude.  There  is,  however,  a  significant  re¬ 
apportionment  of  energy  between  the  frozen-flow  and  thermal 
loss  categories.  This  is  largely  because  of  the  increase  in  m  of 
the  higher  power  thruster  and  the  corresponding  decrease  in 
flow  transit  time  relative  to  molecular  transport  times.  Similar 
trends  are  anticipated  for  other  more  complex  propellants  such 
as  N2H4.  The  loss  categories  with  the  greatest  potential  for 
improved  performance  are  frozen-flow  and  thermal.  Frozen- 
flow  losses  may  be  further  split  into  energies  associated  with 
translation,  rotation,  vibration,  electronic  excitation  (consid¬ 
ered  negligible),  dissociation,  and  ionization.  The  thermal  cat¬ 
egory  includes  radiative  losses  from  the  plasma  and  the  anode 
as  well  as  conduction  to  the  anode  and  supporting  structures. 
Finally,  the  thrust  category  is  taken  to  be  only  the  directed 
kinetic  energy  of  the  arcjet.  If  velocity  profile  losses  are  small, 
its  value  should  be  similar  in  magnitude  to  the  conventional 
thrust  efficiency. 


■  Translation 
Q  Rotation  liVIbratlon 
□  loniratlon 
M  Dissociation 


Frozen  Row 


Thermal 


B  Thruet 


Fig.  12  Output  energy  distribution  for  low  and  intermediat 
power  hydrogen  arcjets. 


Achieving  control  of  the  mass  and  energy  transport  betweei 
the  plasma  and  the  anode  is  a  key  element  of  this  strateg) 
The  difficulty  that  arises  is  that  the  physical  processes  contrib 
uting  to  the  loss  in  performance  are  strongly  coupled.  For  ex 
ample,  the  plasma  residence  time  can  be  increased  in  high 
density  regions  of  the  arcjet  with  a  dramatic  increase  in  ioni 
and  atomic  recombination  rates  (zone  3,  Fig.  1).  Howevei 
without  special  consideration  given  to  the  anode  design,  littl 
of  the  additional  chemical  energy  released  will  be  convene 
to  useful  thrust.  Instead,  it  will  be  transported  by  molecula 
diffusion  to  the  anode  surface  and  realized  as  thermal  losj 
Increased  thermal  loads  resulting  from  this  process  reduc 
thruster  lifetime. 

The  ability  to  manage  the  plasma-anode  energy  transfer  re 
quires  an  understanding  of  arcjet  plasma  physics,  including  th 
dynamics  of  the  arc  attachment.  Unfortunately,  the  hostil 
physical  environment  and  small  size  limit  attempts  to  measur 
the  state  of  the  interior  flow.  The  bulk  of  the  diagnostics  t 
date  involve  measurements  of  external  flowfield  parameter 
chemical  composition,  and  surfaces  properties.  Internal  d> 
namics  have  been  left  to  numerical  models.  The  logic  for  thi 
approach  rapidly  becomes  circular  since  the  confidence  in  sue 
predictions  rests  on  the  fidelity  of  the  measurements  used  t 
validate  them.  Much  remains  to  be  done  to  develop  soun 
engineering  design  tools  to  transition  from  concept  to  practice 

Approach  to  Performance  Model  Development 

The  use  of  sophisticated  numerical  techniques  to  predic 
arcjet  performance  has  grown  considerably  with  the  transitio 
of  the  arcjet  into  flight  application.  Early  work  employing  sirr 
plified  plasma  dynamics  and  equilibrium  chemistry*^'®  ha 
been  replaced  by  full  Navier-Stokes  solvers.  These  ap 
proaches  typically  include  finite  rate  chemistry  and  transpo: 
prop)erties.  Maxwell’s  equations,  and  some  model  for  energ 
transport  across  electrode  boundaries.  At  the  present  time,  sir 
gle-fluid,  single-temperature  (1-F)  models  represent  state-oi 
the-art  design  capability.*’  **  Such  models  have  shown  reasor 
able  agreement  with  exp>eriment,  particularly  when  th 
plasmadynamic  solution  is  coupled  with  the  anode  thermal  di‘ 
tribution.  The  self-consistent  treatment  of  the  anode  an 
plasma  is  essential  at  high  Pn^ni  where  material  failure  is  ar 
ticipated.  Any  numerical  evaluation  of  novel  high-performanc 
concepts  without  this  capability  must  be  held  suspect. 

Figures  13  and  14  show  a  comparison  of  a  1-F  model 
predictions  with  measured  values  for  velocity  and  tempjeratui 
of  excited  state  H  atoms  at  the  exit  plane  of  a  low-power  T 
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Fig.  13  Comparison  of  measured  and  predicted  exit  plane  veloc¬ 
ity  profiles. 
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Fig.  14  Comparison  of  measured  and  predicted  exit  plane  tem¬ 
perature  profiles. 

arcjet.*®  The  arcjet  operating  conditions  were  Pn  =  1.5  kW  and 
m  =  13.1  mg/s.  The  agreement  between  predictions  and  mea¬ 
surements  for  both  velocity  and  temperature  is  remarkably 
good.  Similar  agreement  has  been  demonstrated  for  P tilth 
ranging  from  50  to  125  MJ/kg  (i.e.,  a  factor  of  2.5  in  Pn). 
Although  such  results  are  encouraging,  the  fully  coupled  1-F 
approach  is  not  without  its  caveats.  The  validity  of  transport 
properties  and  chemical  rate  mechanisms  at  the  plasma  tem¬ 
peratures  are  only  the  most  obvious  of  the  uncertainties.  Per¬ 
haps  the  most  serious  shortcoming  of  the  1-F  models  is  that 
they  do  not  accurately  represent  the  electron  transport  or  chem¬ 
ical  activity  in  cooler  regions  of  the  plasma.  This  results  in 
low  estimates  for  electron  population  densities  and  cr,  and 
overestimates  the  ohmic  dissipation  near  the  anode.  To  com¬ 
pensate,  some  functional  dependence  of  cr,  on  Pn  and  other 
plasma  parameters  must  be  assumed  in  the  near-anode  regions. 
The  drawback  of  this  approach  is  that  the  resulting  perfor¬ 
mance  predictions  may  only  be  valid  for  a  limited  range  of 
electrode  geometries  and  arcjet  operating  conditions. 

In  principle,  a  more  general  and  physically  correct  repre¬ 
sentation  of  the  plasmadynamics  can  be  obtained  using  mul¬ 
tifluid  models.^®”^^  In  this  approach,  the  different  species  are 
generally  grouped  into  particle  classes  with  similar  tempera¬ 
tures  (energies)  and  velocities.'^  The  dynamics  of  each  particle 
class  are  then  considered  to  evolve  indef>endently  with  the 
exchange  of  mass,  momentum,  and  energy  between  classes 
expressed  by  appropriate  models. 

The  choice  of  class  structure  depends  on  the  nature  of  the 
problem  and  the  detail  to  which  the  solution  must  be  known. 
For  example,  1-F  models  provide  sufficiently  accurate  perfor¬ 
mance  trends,  but  do  not  properly  represent  the  cooler  regions 
of  the  plasma.  A  two-fluid,  two-temperature  model  with  heavy 


particle  and  electron  classes  will  better  represent  the  cooler 
regions,  but  cannot  capture  the  physics  of  the  electrode  pre¬ 
sheath  or  primary  sheath  regions.  Consideration  of  the  ion  and 
electron  dynamics  at  the  plasma- anode  interface  is  required 
to  reach  this  level  of  detail.  To  accurately  estimate  the  full 
thermal  loading  of  the  anode,  therefore,  a  multi  fluid  perfor¬ 
mance  model  must  include  neutral,  ion,  and  electron  classes. 

The  more  powerful  multi  fluid  approach  provides  a  clearer 
view  of  the  important  physical  processes  and  considerable 
modeling  flexibility.  However,  it  brings  with  it  many  of  the 
difficulties  of  the  1-F  models  in  addition  to  several  others.  The 
determination  of  appropriate  transport  properties  and  chemical 
rates  remain  important  issues.  In  the  multifluid  case,  however, 
species  production  and  destruction  must  also  be  considered  in 
terms  of  the  exchange  of  mass,  momentum,  and  energy  be¬ 
tween  classes.  Although  the  destination  class  of  the  product 
species  is  clear,  the  changes  in  class  momentum  and  energy 
depend  on  the  choice  of  reaction  model. 

The  possibility  of  mass  exchange  between  classes  leads  to 
one  of  the  difficulties  in  applying  this  method,  cases  that  in¬ 
volve  a  vanishing  class  population.  Consider  the  situation  at 
the  root  of  the  arc  where  the  plasma  is  nearly  fully  ionized  or 
along  the  anode  wall  near  the  exit  plane  where  the  ion  and 
electron  populations  are  considerably  diminished.  Under  these 
conditions,  the  intraclass  collision  frequencies  become  so  low 
that  the  classwise  transport  of  mass  and  energy  can  no  longer 
be  based  on  continuum  assumptions,  for  example,  the  Na- 
vier- Stokes  equations.  Note,  however,  that  the  minority  class 
properties  naturally  limit  to  those  of  the  majority  class  as  a 
result  of  intraclass  exchange  processes.  Multifluid  solutions  re¬ 
quire  approximations  for  the  minority  class  dynamics  that  limit 
in  a  physically  consistent  way  as  their  population  becomes 
small.  A  minimum  number  density  or  floor  can  be  imposed 
throughout  the  plasma,  but  a  value  must  be  chosen  a  priori 
that  is  small  enough  to  represent  the  physics,  yet  large  enough 
to  avoid  physically  implausible  results.'® 

Finally,  the  description  of  mass  and  energy  transport  at  the 
electrode  boundaries  requires  appropriate  boundary  conditions 
for  each  of  the  particle  classes.  The  motion  of  charged  particles 
remains  coupled  to  neutrals  to  the  extent  that  the  plasma  re¬ 
mains  collisional.  If  the  electrode  surface  is  chemically  active, 
the  surface  recombination  rate  of  ions  and  electrons  also  sets 
a  boundary  condition  on  the  production  of  neutral  class  spe¬ 
cies.  The  choice  of  surface  chemistry  model  will  therefore 
have  a  strong  influence  on  estimates  for  current  attachment  and 
energy  transport  across  the  electrode  surface.^^  In  practice, 
such  boundary  conditions  have  proven  to  be  stiff  and  difficult 
to  implement.  Although  the  use  of  the  number  density  floor 
mentioned  earlier  may  alleviate  this  problem,  a  perhaps  better 
alternative  would  be  to  develop  integral  models  of  the  near¬ 
electrode  region. 

Considering  these  modeling  options,  a  development  path 
that  builds  upon  the  strengths  of  both  the  single-  and  multifluid 
models  may  be  advanced.  The  general  objective  would  be  to 
develop  an  understanding  of  internal  arcjet  dynamics  sophis¬ 
ticated  enough  to  develop  and  evaluate  competing  design  op¬ 
tions.  Experimentally  validated  1-F  approaches  would  continue 
to  provide  guidance  for  experiments  and  the  evaluation  of  new 
electrode  configurations.  Their  application  would  provide 
timely  and  accurate  performance  predictions  for  propellants 
and  operating  conditions  appropriate  for  advanced  hold,  lift, 
and  move  missions.  This  will  require,  at  a  minimum,  a  much 
improved  description  of  the  plasma  physics  in  regions  near  the 
electrodes. 

Multifluid  modeling  techniques  should  be  developed  to  de¬ 
scribe  the  arcjet  plasmadynamics  from  a  first  principles  point 
of  view.  These  detailed  multifluid  models,  along  with  comple¬ 
mentary  diagnostics,  will  serve  as  investigative  tools  to  iden¬ 
tify  and  to  demonstrate  an  understanding  of  the  critical  per¬ 
formance  issues.  An  obvious  near-term  goal  of  this  effort 
would  be  to  develop  physical  models  that  relate  the  far-field 
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plasma  conditions,  current  density,  and  surface  reactivity  to  the 
anode  thermal  loading.  This  work  would  contribute  to  the  de¬ 
velopment  of  a  generic  wall  model  to  improve  single-fluid  per¬ 
formance  estimates. 

Approach  to  Arcjet  Diagnostic  Development 

The  development  of  high-fidelity  diagnostics  to  interrogate 
both  the  arcjet  plasma  and  electrode  surfaces  is  critical  to  the 
validation  of  numerical  models  as  well  as  the  measurement  of 
advanced  concept  performance.  The  basic  suite  of  measure¬ 
ments  taken  during  an  arcjet  life  test  include  current,  voltage, 
propellant  feed  pressure,  mass  flow  rate,  temperature,  and 
thrust.  From  these  parameters  one  then  derives  /^p  and  17  and 
a  qualitative  measure  of  thruster  health.  It  is  clear  that  such 
data  do  not  provide  detailed  information  concerning  energy 
and  mass  transport  within  the  arcjet.  Recent  efforts  have  been 
made  to  employ  advanced  diagnostic  techniques  with  experi¬ 
ment  setups  carefully  designed  to  isolate  the  critical  phenom¬ 
ena.  Examples  include  the  manufacture  of  anodes  especially 
configured  to  measure  current  distributions,^'*  internal  plasma 
temperatures,^^  spectroscopic  measurements  of  excited  state 
species  in  the  plume,^^"”  and  particle  fluxes  in  the  plume.^ 
Preferred  measurement  techniques  do  not  alter  the  flowfield 
or  electrical  characteristics  of  the  plasma,  are  independent  of 
electrode  configuration,  and  can  be  used  over  a  wide  range  of 
operating  conditions.  Such  measurements  should  determine  the 
distribution  of  energy  among  the  different  loss  modes  as  a 
function  of  arcjet  operating  condition.  This  provides  a  means 
to  evaluate  global  changes  in  performance  that  result  from  de¬ 
sign  changes.  For  example,  the  simultaneous  measurement  of 
frozen-flow  and  thermal  losses  at  high  Pnim  is  essential  since, 
at  least  for  conventional  geometries,  they  trade  in  an  almost 
one-to-one  fashion.  This  is  a  formidable  objective  since  it  re¬ 
quires  the  measurement  of  species  concentrations,  tempera¬ 
tures  or  energies,  and  velocities  as  a  function  of  operating 
condition. 

A  partial  view  of  the  frozen-flow  and  thermal  loss  dynamics 
can  be  constructed  by  combining  recent  results  of  emission 
spectroscopy,  laser- induced  fluorescence  (LIF),  and  thermal 
imaging  measurements.  Emission  spectroscopy  has  been  used 
to  identify  major  contributors  to  frozen-flow  losses.  For  ex¬ 
ample,  it  can  provide  estimates  of  //*  number  density  at  the 
exit  plane  of  H2  arejets.*^  This  data  can  then  be  related  via 
nonequilibrium  rate  models  to  the  local  electron  and  ion  num¬ 
ber  densities.  For  sufficiently  high  number  densities.  Stark 
broadening  of  the  excited  state  emission  can  also  provide  a 
more  direct  estimate  of  the  electron  or  ion  number  density. 
Having  arrived  at  this  estimate,  some  measure  of  the  exit  plane 
velocity  distribution  is  neeeded  to  determine  the  local  ion  flux 
and  therefore,  the  total  ionization  loss.  The  first  experimental 
measurements  of  low-power  H2  arcjet  exit  plane  velocities 
have  been  obtained  for  H  atoms  over  a  range  of  Pnlrh  using 
LIF  techniques,'^  These  velocity  profiles  can  be  used  with  the 
number  density  data  to  estimate  ionization  losses.  Perhaps 
more  importantly,  the  fact  that  1-F  arcjet  performance  model 
predictions  correlate  reasonably  well  with  these  measurements 
(see  Fig.  13)  allows  the  extension  such  loss  estimates  to  con¬ 
ditions  where  LIF  data  are  sparse  or  nonexistent. 

As  shown  in  Fig.  12,  the  energy  tied  up  in  ionization  rep¬ 
resents  only  a  small  fraction  of  the  total  frozen-flow  loss  in 
hydrogen  arejets.  The  fact  that  the  ion  population  also  in¬ 
creases  as  the  total  frozen-flow  loss  increases  offers  the  hope 
that  it  can  be  taken  as  a  trend  indicator.  For  example,  one  could 
hypothesize  that  an  electrode  configuration  that  produces  lower 
ionization  losses  also  produces  lower  total  frozen-flow  losses. 
This  requires  the  measurement  and  correlation  of  the  major 
species  number  densities  with  model  predictions.  In  particular, 
the  primary  objective  should  be  the  determination  of  the  num¬ 
ber  densities  of  the  neutral  products  of  dissociation,  since  they 
represent  the  largest  contributors  to  the  dissociation  losses.  Ef¬ 
forts  toward  this  objective  have  so  far  met  with  only  limited 
success.^ 


The  measurement  of  the  interior  surface  thermal  loading  an* 
the  resultant  thermal  losses  presents  another  technical  chal 
lenge.  Some  success  has  been  achieved  using  full  thermr 
models  of  the  arcjet.  This  method  typically  assumes  a  certai 
thermal  loss  and  iterates  upon  the  spatial  variation  of  cnerg 
flux  on  interior  surfaces  to  achieve  a  match  between  the  pre 
dieted  and  measured  exterior  surface  temperatures.  When  a 
acceptable  agreement  has  been  obtained,  one  can  estimate  th 
total  conductive,  connective,  and  radiative  energy  exchang 
associated  with  the  exterior  and  interior  surfaces.  The  shor 
comings  of  this  approach  are  twofold.  First,  thermocouple 
cannot  survive  the  arejet’s  high  steady-state  anode  surface 
temperatures.  The  lack  of  temperature  data  greatly  increase 
the  uncertainty  of  the  loading  estimates  since  radiation  froi 
the  anode  exterior  is  the  largest  contributor  to  the  thermal  lo^ 
category.  Second,  and  most  important,  the  exterior  surfac 
temperature  profile  is  not  uniquely  determined  by  the  interic  | 
surface  fluxes,  i.e.,  the  exterior  temperatures  are  insensitive  1 
variations  in  the  assumed  internal  loading  profiles. 

The  first  shortcoming  has  been  addressed  in  most  respec 
by  thermal  imaging  of  the  anode  exterior  surfaces."  '^  Thee 
techniques  provide  a  nearly  continuous  measure  of  the  surfac 
temperatures  and  greatly  increase  the  sensitivity  of  the  iterati\ 
scheme  described  earlier.  If  necessary,  further  refinements  i 
the  loading  distribution  may  be  obtained  using  profiles  pn 
dieted  by  numerical  performance  models. 

The  second  shortcoming  can  only  be  solved  by  direct  me: 
surement  of  the  interior  surface  temperature  or  energy  fit 
distribution.  The  development  of  diagnostic  techniques  to  a( 
complish  this  task  should  receive  a  high  priority.  If  arcjet  d( 
signs  are  to  move  toward  significantly  higher  Pnfrh,  the  pn 
cesses  that  occur  at  the  plasma-electrode  boundary  must  t 
better  understood  and  predicted  with  a  reasonable  degree  ( 
confidence.  It  follows  then,  that  nonintrusive  measuremen 
must  be  made  in  the  interior  of  the  arcjet,  or  in  arejet-simil 
plasmas,  to  gain  insight  into  the  critical  physics  and  valida 
models  of  the  key  processes. 

Summary 

Arcjet  technology  offers  significant  benefits  over  other  ele 
trie  propulsion  options  for  many  advanced  mission  applic 
tions.  The  impact  of  the  arcjet  on  hold,  move,  and  lift  missic 
capabilities  will  be  significant  if  next-generation  arejets  ci 
demonstrate  sufficiently  high  A  number  of  design  optioi 
were  discussed,  including  operation  at  decreased  mass  flu 
increased  power,  and  increased  efficiency.  It  appears  that  tl 
only  realistic  way  to  achieve  the  desired  /,p  is  to  improve  upc 
present  thrust  efficiencies  by  reducing  both  frozen-flow  ai 
thermal  losses.  The  objective  is  to  reduce  surface  energy  de 
sities  by  controlling  the  transport  of  mass  and  energy  to  tl 
anode  surface.  Improvements  should  include  the  continued  r 
finement  and  validation  of  single-fluid  models  to  provii 
timely  performance  estimates,  the  use  of  first  principles  mi 
tifluid  models  to  investigate  critical  transport  processes,  tl 
extension  of  these  models  to  general  performance  estimate 
and  the  development  of  diagnostic  techniques  to  query  t; 
plasma  and  interior  surfaces.  In  this  endeavor,  the  approach 
to  measurement  and  modeling  development  are  intimately  co 
pled. 

Adopting  this  approach  has  the  added  benefit  that  the  arej 
TIP  and  lifetime  advantage  over  other  electric  propulsion  d 
vices  remains  unchallenged.  Moreover,  improvements  in  bas 
process  efficiencies  will  build  directly  on  efficiency  gain 
from  improvements  in  electrode  material  properties  and  a 
vanced  schemes  to  process  high-energy  fluxes.  The  desi; 
principles  derived  from  this  effort  will  apply  to  all  classes 
arejets  and  all  propellants. 
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